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ABSTRACT 


The  purpose  of  this  work  was  to  identify  and  study  pulse 
compression  techniques  useful  in  LPI  (l.e.,  Low  Proabability  of  J 

Intercept)  applications.  More  specifically,  the  technique  pro- 
posed  should  be  capable  of  compression  ratios  of  from  128  to 
1024  or  higher,  and  Integrated  si delobes  of  -35  to  -30  db 
throughout  the  compression  range;  furthermore,  it  should  be  easy 
to  mechanize,  and  should  result  in  minimal  S/N  loss.  It  was 
concluded  that  nonlinear  FM  is  the  pulse -compression  scheme 
capable  of  satisfying  all  these  requirements,  and  existing 
technologies  make  its  mechanization  feasible. 
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I.  INTRODUCTION 

The  probability  of  detection  of  a radar  signal  in  the  presence  of  white 
Gaussian  noise  is  maximized  if  the  receiver  is  a matched  filter  to  the  signal 
[1].  In  this  case  the  receiver  response  becomes  the  signal  autocorrelation 
function.  In  certain  radar  applications.  It  Is  also  necessary  that  the 
transmitted  waveform  be  rectangular  and  have  as  low  an  amplitude  as  possible 
but  always  consistent  with  the  requirements  for  range  resolution  and  prob- 
ability of  detection.  These  signal  design  constraints  are  best  satisfied  by 
using  pulse  compression  techniques  which  will  result  in  the  minimal  possible 
signal-to-noise  power  loss. 

Nonlinear  FM  signals  have  relatively  low  time  sidelobes  and  result  in  no  loss 
in  signal-to-noisc  ratio  for  as  long  as  the  receiver  is  matched  to  the  trans- 
mitted nonlinear  FM  waveform.  The  compressed  pulse,  however,  is  highly  sensi- 
tive to  frequency  variations,  and  becomes  extremely  distorted  and  develops  high 
sidelobes  as  the  doppler  frequency  is  Increased. 

Since  nonlinear  FM  suffers  no  signal-to-noise  ratio  loss  under  perfect  match 
conditions,  it  becomes  a serious  candidate  for  use  In  a pulse  compression  sys- 
tem. In  the  sequel,  a digital  pulse-compression  system  is  proposed  based  on 
nonlinear  FM  and  capable  of  operating  over  an  arbitrarily  wide  band  of  doppler 
frequencies  without  violating  the  requirements  of  specified  sidelobe  level  and 
signal-to-noise  ratio  degradation. 


11.  NONLINEAR  FM  PULSE  COMPRESSION  SIGNALS 

Consider  a narrow  band  signal  s(t)  centered  at  a carrier  frequency  f£  with 
an  envelope  a(t),  and  a phase  modulation  given  by  t(t).  That  Is, 


s(t)  * a(t)cos[2irf  t + <>(t)3  • -T/2  < t < T/2 

V * "■ 

* 0,  elsewhere  * (1) 

• . 

where  T Is  the  pulse  length.  The  equivalent  complex  representation  of  the 
modulation  Is  given  by 


a(t)  * a(t)  exp[j$(t)]  (2) 

and  its  corresponding  Fourier  transform  Is  expressed  as 

A(f)  * a(t)  exp[j(4>(t)  - 2irft}]dt 

* A(f)  exp[je(f)]  (3) 

where  A(f)  denotes  the  modulus  of  A(f).  In  a similar  manner,  one  may  also 
write 

a(t)  = j A(f)  exp[j(e(f)  + 2nft}]df 

c a(t)  exp[j<j>(t)]  •"  (3a) 


where  a (t ) is  the  modulus  of  a(t). 

The  moduli  a(t)  and  A(f)  are  related  via  Parseval's  theorem  which  Is 
analytically  expressed  by 

J"  a2(t)dt  = jV(f)df  = 2E 

where  E is  the  energy  in  the  real  signal. 


(4) 


The  signal  autocorrelation  function,  z(t),  which  appears  at  the  output  of  the 
matched  filter  receiver,  is  given  by 

z(t)  a j a(i)  a*(t+T)  dx  (5) 
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where  a*(t)  Is  the  complex  conjugate  of  a(t).  Using  the  Fourier  transform 


of  a(x),  one  may  rewrite  Eq.  (5)  as 

z(t)  » J A2(f)  exp[j2irft]df 

(6) 

or 

A(f)  *|  fz(t)  exp[-j2*ft]dtj  ^ 

m 

The  general  problem  which  includes  the  nonlinear  FM  case  as  a special  case 
may  be  stated  as  follows: 

"Given  the  two  moduli  a(t)  and  A(f),  is  it  possible  to  obtain  a pair  of  phase 
functions  $(t)  and  e(f)  such  that  equation  (3)  (or  3a)  is  satisfied?" 


With  the  exception  of  very  few  well-known  cases,  the  general  exact  solution 
to  this  problem  for  arbitrary  functions  a(t)  and  A(f)  is  not  known.  This 

stems  from  the  fact  that  Fourier  theory  does  not  provide  a general  set  of 

necessary  and  sufficient  conditions  which  must  be  satisfied  by  the  moduli  of 
a potential  Fourier  transform  pair,  that  is,  by  a(t)  and  A(f)  In  this  case. 

(Note,  once  again,  that  the  reference  here  is  to  a(t)  and  A(f)  and  not  to 
£(t)  and  Nevertheless,  one  may  still  insist  on  an  approximate  solu- 
tion and  accept  it  or  reject  it  depending  on  how  well  it  satisfies  expressions 

(3)  and  (3a).  Such  approximate  solutions  may  be  pursued  by  making  use  of  the 

« 

stationary  phase  method  of  approximate  integration  [2].  This  method  may  be 
applied  either  to  expression  (3)  or  (3a). 


According  to  the  principle  of  stationary  phase,  the  major  contribution  to  the 
Integral  of  a rapidly  oscillating  function  will  result  from  regions  over  which 
the  phase  of  the  integrand  functions  is  "stationary";  that  is,  from  regions 
over  which  the  derivative  of  the  phase  function  is  zero. 
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Applying  this  technique  to  expression  (3),  one  has 
^ [*(t)  - 2 it  ft]  = 0 

or  ♦,(t)  * 2*f  (8) 

Expression  (8)  implies  that  the  points  in  time,  t,  at  which  the  total  phase 
becomes  "stationary"  will  depend  on  the  frequency  f. 

(toe  obtains  a similar  relationship  using  expression  (3a).  That  is, 

[0(f)  + 2irft]  =0 

or  6* (f)  * - 2wt  (9) 

Expression  (9)  implies  that  the  frequencies,  f,  at  which  the  total  phase 
becomes  "stationary"  will  depend  on  the  time,  t. 

It  becomes  apparent  that  in  order  to  obtain  $(t),  it  is  necessary  that  we 
determine  the  functional  dependence  of  f on  t,  or 

f a g(t)  (9a) 

Use  of  (9a)  in  (8)  and  subsequent  integration  will  provide  the  phase  $(t). 

Similarly,  one  may  use  the  inverse  of  (9a),  that  is, 

t • g'^f)  (9b) 

and  expression  (9)  to  determine  e(f). 

Determining  the  functional  dependence  expressed  by  (9a)  or  (9b)  is  the  ulti- 
mate objective  of  the  stationary  phase  approximation  method.  This  is  illu- 
strated by  considering  expression  (3a)  and  expanding  its  phase  function  in 
a 3-term  Taylor's  series  about  the  stationary  point,  say,  f = e,and  performing 
the  Indicated  operations.  Key  et  al.  [3]  have  shown  that  this  yields  the 
approximate  expression 

/ 
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a(t) » 


MsL 


exp[j{27ret  + 0(e)  + ir/4}] 


00) 


■ 

where  (+)  Is  associated  with  e"(e)  > 0 arid  (-)  with  e"(e)  < 0.  Expression 
(10)  for  a(t)  may  be  separated  into  the  modulus  and  phase  functions  as  follows: 


a(t)  = /hr 


and 


V |eV)| 

$(t)  = 2irft  + 0(f)  + tt/4 


(ID 

(12) 


where  e,  denoting  the  stationary  point  in  frequency,  has  been  replaced  by  f 
for  clarity  and  in  being  consistent  with  expressions  (8)  and  (9). 


Differentiation  of  expression  (9)  and  subsequent  substitution  in  (11)  yields 
a2(t)  dt  = A2(f)  df  (13) 

which  in  a sense  establishes  in  differntial  form  the  dependence  of  f on  t or 
vice  versa.  Furthermore,  expression  (13)  prescribes  the  manner  in  which  the 
moduli  a(t)  and  A(f)  enter  in  this  dependence. 


According  to  Parseval’s  theorem  (see  expression  (4)),  Eq.  (13)  must  yield 


J a2(t)  dt  = | A2(f)  df 

(14) 

But  rore  importantly,  given  time  t (or  frequency  f)  we  want 
sponding  frequency  f (or  time  t)  such  that 

to  know  the  corre- 

j 

1 

* 

9 

| 

1 
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< 

J a2(y)dy  = j AZ(x)dx 

(15) 

- 

or  equally  valid 

t _ 

* 

aZ(y)dy  = AZ(x)dx 

(16) 
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Equations  (15)  and  (16)  provide  two  possible  solutions  regarding  the  de- 
pendence of  f on  t,  or  vice  versa,  given  two  moduli  functionsa(y)  and  A(x). 

If  the  integrations  implied  by  (15)  and  (16)  are  not  tractable  analytically, 
one  may  employ  computational  techniques  to  obtain  the  dependence  of  t on  f 
or  vice  versa. 

One  special  case  of  interest  is  the  one  in  which  the  signal  envelope  or 
modulus  a(t)  is  given  by 

a(t)  • 1,  -T/2<t<.T/2 
= 0,  elsewhere 

and  the  spectral  function  or  modulus  A(f)  is  such  that  when  processed  by  the 
matched  filter  receiver  results  In  the  prescribed  response.  The  choice  of 
A(f)  for  a specified  si  delobe  level  and  compressed  pulse  shape  can  be  made 
from  any  of  several  well-known  families  of  spectral  functions  [3,4]. 

Appendix-1  provides  the  computer  program  used  to  compute  numerically  the  re- 

quired  phase  function  q>(t)  given  a (x)  = 1,  -T/2  <_  t <_T/2,  and  arbitrary 

o 

spectral  functions  for  A (f).  In  the  program,  use  of  expression  (15)  is  made 
to  compute  the  value  of  f corresponding  to  a given  t,  and  this  information 
is  subsequently  utilized  in  Eq.  (8)  to  evaluate  the  derivative  <|>'(t)  at  that 
point  in  time.  From  the  numerically  established  behavior  of  the  derivative, 
the  program  proceeds  and  evaluates  <j>(t)  which  is  the  phase  to  be  assigned  to 
the  transmitted  signal. 


--  H 
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III.  NONLINEAR  FM  SIGNALS  IN  THE  PRESENCE  OF  DOPPLER  SHIFT 
The  compressed  pulse  resulting  from  nonlinear  FM  pulre-compression  systems 
is  highly  sensitive  to  frequency  shifts  due  to  doppler.  This  is  due  to  a 
resulting  phase  mismatch  condition  which  is  responsible  for  causing  severe 
distortions  in  the  receiver  response  both  in  terms  of  degradation  in  the 
main  lobe  and  the  level  of  the  range  si  delobes. 

The  received  spectrum  is  given  by 

S(w)  = A (w+uj)  expCjetu+uj)]  (17) 

where  ■ 2*^.  the  doppler  frequency  shift  in  rad/sec.  The  matched-filter 
transfer  function  is  given  by 

H(u)  = A(oj)  exp[-je(u)]  • (18) 

and  the  output  of  the  matched  filter  by 

G(u)  = ACu+Wj)  A(w)  expCjOCu+u^i-atu) }]  (19) 

The  function  e (to)  is  obtained  using  procedures  developed  in  the  preceding 
section.  _ ~ 


In  order  to  overcome  the  detrimental  effects  of  phase  mismatch  due  to  doppler 

1 

frequency  shifts,  it  is  proposed  that  a bank  of  filters  be  designed  which  are 
perfectly  matched  in  phase  to  signals  received  at  equally  spaced  frequencies 


spanning  the  entire  Interval  of  doppler  frequencies  of  interest  (i.e. , 

Vfd  to  W- 

C UMaX  C uM«x 


The  center  frequencies  of  the  matched  filters  used  in  adjacent  doppler 
channels  are  so  spaced  that  the  resulting  deterioration  in  range  sidelobe 


levels  is  kept  within  prespecified  bounds. 


The  design  of  matched  filters  for  processing  nonlinear  FM  signals  is  a very 
difficult  technical  undertaking,  and  the  task  becomes  even  more  complex  if 


-7- 


the  filters  are  required  to  process  a wide  variety  of  spectral  functions, 
A(f).  Such  a need  may  arise  either  from  technical  or  from  tactical  con- 
siderations, and  it  is  best  satisfied  by  implementing  the  matched  filters 
and  associated  signal  processing  numerically.  Such  an  approach  requires 
that  we  design  the  matched-filter  transfer  function  in  a way  which  accounts 
for  phase  distortions  due  to  doppler  frequency  shifts.  Knowledge  of  the 
dependence  of  A(to)  on  u>  facilitates  accounting  of  the  doppler  effects  on 
the  signal  and  its  subsequent  processing.  In  order  to  obtain  the  filter 
characteristics  for  perfect  match  in  the  presence  of  a doppler  frequency 
shift  equal  to  one  must  obtain  e(u>)  from  the  moduli  a(t)  and  AU+u^). 
The  modulus  a(t)  being  the  envelope  of  the  transmitted  signal  is  assumed 
to  remain  fixed,  but  A(w)  is  adjusted  to  reflect  the  frequency  shift  due 
to  doppler. 

The  required  spacing  between  the  center  frequencies  of  adjacent  doppler 
filters  for  a specified  sidelobe  or  sianal-to-noise  ratio  deterioration 
is  determined  via  numerical  techniques  utilizing  the  program  given  in 
Appendix-1.  In  the  case  of  a square-root  Taylor  frequency  spectrum, 
adjacent  doppler  channels  with  their  center  frequencies  placed  at  points 
corresponding  to  a target  velocity  of  1.76  of  MACH-1  result  in  a signal- 
to-noise  ratio  deterioration  of  less  than  -0.22  db;  the  corresponding 
distortion  in  the  time  response  yields  -18.1  db  sidelobe  level  instead 
of  the  design  value  of  -33.2  db. 


Figure  1 offers  a processing  scheme  for  non-linear  FM  signal  waveforms. 


Nonlinear  FM  Pulse  Compression  System 


IV.  DESIGN  CONSIDERATIONS  FOR  THE  NON-LINEAR  FM  PULSE  COMPRESSION  SYSTEM 
In  order  to  accommodate  information  rates  and  signal  waveforms  compatible 
with  modern  radar  system  requirements,  the  design  of  the  nonlinear  FM 
pulse  compression  system  must  display 

a/  flexibility  In  the  choice  of  transfer  function,  and 
b/  bandwidth  capability  of  the  order  of  many  megahertzs. 

The  issue  of  system  flexibility  is  easily  settled  by  assuming  digital  imple- 
mentation of  the  processor.  This  choice  usually  puts  extra  pressure  on  the 
system  designer  to  develop  FFT  (Fast-Fourier  Transform)  units  capable  of 
coping  with  large  bandwidths.  In  the  entire  system  shown  in  Fig.  1,  the  FFT 
units  are  the  only  design  challenges,  and  with  present-day  technology  and 
existing  FFT  algorithms  [5-12]  their  implementation  has  become  both  economically 
and  technically  feasible.  The  factors  which  determine  the  size  of  the  FFT 
unit,  and  thus  its  cost,  are  the  time-bandwidth  product,  speed  and  precision 
requirements  of  the  system. 

At  this  point  in  time,  the  design  of  FFT  processors  has  reached  sufficient 
maturity  to  satisfy  most  of  the  signal  processing  needs  encountered  in 
modern  radar  systems.  The  analysis  and  error  sensitivity  of  various  FFT 
algorithms  have  been  well  documented  [13-20],  and  there  Is  no  need  for  their 
duplication  here.  It  should  be  pointed  out,  however,  that  new  advances  In 
the  area  -of  parallel  computational  structures  [21-25]  point  to  more  cost- 
effective  FFT  processor  designs  in  the  future,  and  thus  to  more  economical 
and  flexible  signal  processing  systems. 
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APPENDIX-1 


The  program  given  In  this  Appendix  uses  square-root  Taylor  weighting  In 
the  frequency  domain  for  -40  db  sidelobe  level.  The  Input  quantities  DELI 
and  DELTA  ( ) specify  the  doppler  frequency  shift  we  wish  to  center  the 
matched  filter  under  design  at,  and  the  sequence  of  doppler  frequencies  for 
use  In  testing  the  performance  of  the  matched  filter,  respectively.  Note 
that  DELTA  or  DELI  are  normalized  quantities  corresponding  to  the  parameter, 

6 = (2  u/c)(fc/Af) 

where  u * target  radial  velocity;  c = velocity  of  light;  f * carrier 
frequency;  and  Af  = signal  bandwidth.  Furthermore,  many  of  the  parameters 
are  assigned  values  within  the  program  and  not  thru  a READ  statement. 

ND  * Number  of  doppler  frequencies  to  be  investigated 
N * Number  of  sample  points  to  be  processed  j M = n/i6. 

NSTAGE  = A positive  integer  such  that  N = 2 **  (NSTAGE).  It  denotes  the 
nunber  of  FFT  stages. 

AF ( ) * Coefficients  determining  the  polynomial  signal  taper  employed.  The 
illustrative  values  used  In  the  program  represent  a -40  db  sii delobe 
square-root  Taylor  weighting.  AF ( ) Is  DIMENSIONED  to  accommodate 
up  to  ten  polynomial  coefficients  (see  p.  20-18  of  RADAR  HANDBOOK 
by  Skolnlk.) 

D * Compression  ratio.  The  compression  ratio  is  set  equal  to  (N/4  - 0.25) 
which  Indicates  four  samples  per  cycle  with  the  first  sample  at 
time  t ■ 0,  and  the  last  sample  omitted.  (Also,  D ■ C0MPR) 

MF  ■ N/4 


I 
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MPSI  * Number  of  Fourier  coefficients  wanted  for  approximating  the  signal 
phase  function  PSI ( ) in  the  time  domain. 


DELI  = Center-frequency  at  which  the  nonlinear  filter  is  tuned.  It  is 
given  in  normalized  form  and  has  been  defined  earlier  in  the 
Appendix. 

DELTA(  ) = Doppler  frequencies  expressed  in  normalized  form  and  in  a manner 

I similar  to  DELI . 

I 

i 

Discussion  of  subroutines  and  functions 

1.  SUBROUTINE  EVLPSI  (NSTAGE,  N,  NN,  D).  It  evaluates  the  Fourier  co- 
efficients B(  ) for  the  approximation  of  the  real  signal  phase  function 
PSI(t). 

\ 2.  SUBROUTINE  EVFLTR  (NSTAGE,  N,  NN,  D,  FACTOR,  DELI).  It  evaluates  the 

] transfer  function  of  the  nonlinear  filter  tuned  at  DELI.  The  parameter 

FACTOR  is  only  defined  and  computed  within  the  subroutine  for  latter 
use  by  other  subroutines. 


3.  SUBROUTINE  RESPNS  (NSTAGE,  N,  NN,  D,  DEL,  FACTOR, XMAX).  It  evaluates 
the  response  of  the  nonlinear  filter  due  to  a unit-amplitude  rectangular 
pulse  of  real  phase  PSI(t). 

4.  SUBROUTINE  TSOPLT  (MM).  It  provides  a quick  "line-printer"  plot  of  the 
filter  response  (or  of  any  function  computed  within  the  program). 

5.  SUBROUTINE  SAMPLE  (NF,  SIGN,  D,  DEL,  NSTAGE).  It  samples  the  signal 

exp(JPSI(t)) 
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6.  SUBROUTINE  FFT(NSTAGE,  SIGN).  It  evaluates  the  discrete  Fourier 


transform  of  a set  of  sample  points. 

7.  SUBROUTINE  TMPROC  (NF,  0,  DEL,  NSTAGE).  It  is  used  in  conjunction  with 
SUBROUTINE  RESPNS(  ) 

8.  SUBROUTINE  SELECT  (NF,  VARBLE,  OUTPUT).  It  selects  between  the 
functions  SGNL(  ) and  TAPER(  ) depending  on  the  value  NF. 

9.  SUBROUTINE.  INVERT  (NF,  XTOL,  FTOL,  NTOL,  Y).  It  determines  the  root(s) 
of  a function  using  Newton's  method. 

10.  SUBROUTINE  AMPLE  (THIGH,  N,  NN).  It  samples  the  derivative  function  of 
PSI(t). 


11.  SUBROUTINE  INTGRN(NF,  C,  E,  Nl,  NMAX,  STOL,  S).  It  integrates  either 

\ 

function  SGNL  or  TAPER  depending  on  the  value  of  NF. 

12.  FUNCTION  PSI  (D,  DEL,  SIGN,  T).  It  provides  the  real  phase  of  the 
rectangularly  shaped  signal  assuming  square-root  Taylor  weighting. 

13.  FUNCTION  TAPER  (PHI).  It  provides  the  spectral  function  of  the  signal 
which  will  result  in  the  desired  compressed  pulse.  In  this  case,  use 
is  made  of  square-root  Taylor  weighting. 

14.  FUNCTION  SGNL(T).  It  provides  a constant,  unit-amplitude  taper. 

I 

’ 

*- 
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T 


50PLNEW  JOB  DU. D07. AB12Q3, MARINOS, T=  (2,00) ,P=50 
50PLNFW  EXEC  WATFIV,REGIoN»220K  . . , 

SIcoM^E§°J;g5{?^?/MARINOS-TIME“120-PAeES'50 

. DOUBLE  PRECISION  D,DT,DFlOAT 

•*  ?ommHn  pr^t If  ?5at S^SoSef If  a^oU > ,auoi)  »B(ion  >mfour,mf 

COMMON  DATANL(2048) tDATAX(204d) 


DU. 007. AB1 203, MARINOS, T 
W ATF I V , REGI 0N-220K 


=(2,00) ,P=50 


DIMENSION  ‘ , 3?° ^ ®C0N^T^ VIC°NST  jNFLAGj  IFLAG/HjMM  , _ 

MFo&>\x\i\n\  **>'  nstage’  Mj  MPsr~  

103  Cow*  r-  einSSSESE  SBC ~ 

“AF  ( 1 ) = 0 . 3891 154  

AF ( 2 ) =*0 • 0094523 
AF(3)=0. 0048819 

~ ~AF<  4)  =-0.0016105  ' 

AF  ( 5)  =0.0003^74 * ~ 

n=?**nst  age 

D=COVPP  ~ ' ' " " 

PF  ScK?  1.102)  (DELTA  (K)‘,  < = 1 , NDT^  ~ a“i*t '* ^ * 

102  FORMAT  (8F10.0) 

NNs?«N  ~ ' ~ ' 

“*-*•  mfoup=mps  t~~  * — — “ — ' — * — — ' — “*■ — — • — ■ ■-  - * 

Dl=D 

WRITE  (3,1) 

' 1 FORMAT  (I HI)  •’  * - '•  ' 

CALL  FVLPSI  (NSTAGE.N.NNVD) * — 

IF  (IFLAG.GT.O.OR.NFLAG.GT.O)  GO  TO  999 
WRITE  (3.9)  D1 

^9* FORMAT  (3X.31H  THE  COMPR.  RATIONS.  EQUAL  .TO  ITIZ^fFlO. 3,/7)  ; i 
WPITE  (3,302) 

302  FORMAT  ( 3X , • THE  FOURIER  COEFFICIENTS  OF  PSI(T)  ARE',//) 

WRITE  (3,303) 

m format  (6x,sh  n.  • i2«>.8h_b 

DO  304  I=l,MFOUR 
WRITE  (3,305)  1.8(1) 

305  FORMAT  (3X.I4.E20.6./) 

304  CONTINUE  ~ — ~~~ ~ — ~ “rr: — ~~ 

WRITE  (3.1 ) ~ * “ 

CALL  EVFLTR  (NSTAGE. N,NN,0, FACTOR, DELI ) 

_ .INITIALIZE  THE  LOOP  for  the  delta  values.  ....  

15  1FJU.6T.ND)  GO  TO  999 


IF  (J.6T.ND)  GO  TO  999 
DFL=DFLTA(J) 

CALL  PESPNS  (NSTAGE, N,NN,D,DFL*FACT0R,XMAX) 

t.m  T rP  / n a i ft  , '•*  — * • ~ 


' WRITE  (3,9)  D1 
WRITE  (3,21 ) DEL 

21  Format  T 3 X , 2 0 H DELTA  IS  EQUAL  TO 
WRITE  (3.898)  OELl 

)98  FORMAT  < 1 3 X . 3^HF IlTER  CENTERED  A 


.E15.4,//) 


898  FORMAT  < 1 3 X • 2ft  HF 1 LTER  CENTERED  AT  DELTA=  ,E15.4,//) — 

WPITF  (3.98)  XMAX 

98  FORMAT  ( 3X ♦ ’THE  SIGNAL  TO  NOISE  PATIO  IS  EQUAL  TO*,E20.6,«  OBS*. 
1 //) 

WRITE*  (3.99) * 

99  FORMAT  <3X,*THE  AMPLITUDE  OF  THE  RESPONSE  IS  EXPRESSED  IN  DBS1*//) 
CALI.  TSOPLT  (mm) 

WRITE  (3,1)  -**■’ 

GO  TO  15  - 

999  CONTINUE 

eJ8°  COPY  AVIiLIHF  Tn  nnr  now 


COPY  AVAILABLE  TOBDC  DOES  HOT 
PEWIT  FULLY  LEGIBLE  PRODUCTION 


5tre ROUT  I NE'  EVFLTP  ' (NST AGETKl  «NN»  0 * FACTOR  ♦ DEL  1 ) 

COMPLEX  X .DATANL 

DOUBLE  PRECISION  O.DT.PFLOAT 

*1  OOUBLE  PRECISION  AF , XO . GXO ♦ GCONST » WCONST  

COMMON  Dy.X  (?.?04*)  .DATA<20<*6)  .TA(2048>  .A(101>  »B( 
COMMON  OATANL (?049) . D A T A X (2049) 

common  af<io) . xo.gxo.gconsi » wconst » nfl ag. i flag. m, mm 

'COMPUTE  REQUIRED  FILTER  COMPENSATION*  NEXT! 

DT=D/OFLOAT (N-l ) “ 

LST  A6E=NSTAGE* 1 
MF 1 =MF  + 1 

ML.  call  sample  (2««-i..d.o 

call  fft  (LST AGE  t'-l  • ) 

DAT  ANL ( 1 ) =X ( 2 . 1 ) 

' DO  ?50  K=?.MF 

OATANL  (K)=X(2*K)  • 

i "DAT ANL ( MN+2-K ) = X ( 2 .NN* 

?50  CONTINUE 

DO  2S1  K=MF 1 «N 


NSTAGE) 


OATANL (K)=CMPLX (0.0 *0.0)  — — 
— “ DAT  ANL (NN  + 2-K) =CMPLX (O.OtO.O)'  " 
?51  CONTINUE 

DAT ANL (N*l) =CMPLX (0.0. 0.0) 

- —CALL  SAMPLE  1 2.  M • *D«DH.  1 *NST AGE ) 

CALL  F FT  (LSTAGE.-l.) — 

OATANL (1 )=DATANL(1)/X(2.1) 

DC  8001  K=?.MF 

OATANL  (K)  =DATANL  (K)  /X  (2.K)  ' ' 

DAT  ANL (NN  + 2-K) =DATANL (NN*2-K)VX (2 


"CALL  sample  (?.-l.,D.O.*NSTAGEr  

CALL  FFT  (LSTAOEf-1.) 

OATANL (1)=DAT ANL (1)*X(2.1) 

DO  8000  K=2.MF  ~ ~ 

OATANL  (K)  =0ATANL  (i<)  »X  (?.K)  * — 

OATANL  <NN*2-K) =DATANL (NN*2-K> *X (2 .NN*2-K) 

00  CONTINUE 

'THE  COMPLEX  ARRAY  DATANL  CONTAINS  THE  FILTER 


DA  T=  REAL  (OATANL  ( 1 ) ) **2*A  IMAG  (OATANL  ( I ) ) «*2  ' “ — 

DO  8002  K=2.MF 

GaT=  REAL (DATANL (K ) ) **2*A IMAG (OATANL (K) ) **2 

• DAT=DAT*2.*GAT  : "■’* 

02  CONTINUE  • 

THE  CONSTANT  •FACTOR1  IS  EQUAL  TO  THE  DENOMINATOR  IN  THE  DEFINI 
TTON  OF  THE  SIGMAL-TO-NOISE  RATIO. 

--FACTQRsDAT/2. -----  — - 

END1  ' - 


SUBPO'JTINF  QESPNS  (NSTAGE .N.NN.O. DEL .FACTOR. XMAX ) 
COMPLEX  X.DATANL 

DOUBLE  PRECISION  O.DT.DFLOAT”  ' 

PRECISION  AF.XO.GXG.GCONST.WCONST 


DOU°LE 

COMMON 

COMMON 

COMMON 


11 


97 


12 


K=1.NN 

:)=?eal(x(i«k)  ) ««2*aimag(x (i ,k)  j**r 


OT.X (2.2048) .DATA (2040) .TA (2048) »A ( 101 ) . 
OATANL (2048) .DATAX (2048) 

_ AF ( 1 0 ) ♦ XO  *GX0  »GCOMST . WCONST .NFLAG.  IFLAG.M* 
PREPARATION  OF  THE  OUTPUT  OF  THE  RESPONSE 
CALL  TMDROC (P.D.DEL.NSTAGE) 

8atHk?= 

DATA  (K) =DAT A (K) /FACTOR 

TA  (K)  =-n*OFLOAT  (K-l ) *0T 

.CONTINUE  * 

XMAX  = DATA  ( 1 ) — 

DO  97  1=1 .NN 

IF  (DATA(I) .GT.XMAX)  XMAX  = DAT  A ( I ) 

CONTINUE 

DO  1?  K = 1 .MM  - 

DATA  (K) =OATA (N-M*K> 

DATA (K) =10.* ALOGIO  (DATA (K) ) 

TA (K) =TA (N-M*K) 

CONTINUE  — * ~ 

XMAX=10.*ALOG10 (XMAX) 

RETURN 

END  ;r  T ‘ r.-T.TT.^ — vr--?-’:* 


8(101) *MFOUP*MF 

MM” 


TN  FTS7  G N .DTD  ELYUS  TA  G E f 


SU6P0UT  INEZ'S  ample: 

COMPLEX  X.DATANL 

double  precision  d.dt 
‘ “ “ ppeci 

” OT.  _ . 

OAT  ANL ( 20 A8 ) .DaTaX(2048) 

af  ( io)  .xo.gxo.gconst.  ./const,  nflag  .iflag.m.  mm 


GO  TO  100 


3 

00 


DouBl  E 
‘Common 

COMMON 
COMMON 
N=2**NSTAGE 
N1=NM  * 

NM=?*N 

IF  (NF.FQ.l) 

T = -r>/2.-0T 

DO  2 I = J.N  ~”” 

T = T*DT 

Y=PSI (0. DEL. SIGN. T) 
YC=COS(Y)  1 

YS=SlN(Y) 

X(1.I)=CMPLX(YC.YS) 

CONTINUE 

DO  3 I =N1 * NN 

XU  .I)=CMPLX  ( 0.0.0. ~0)' 

CONTINUE 

PFTiicu 

Y = DAT  A ( 1 ) 

YC=COS ( Y ) 

YSsS IN ( Y ) 

X ( 1 . 1 ) =CMRLX ( rC.YS) 

00  4 1=2. MF 
Y=DATA ( I ) 

YC=COS (Y) 

YS=SIN ( Y ) 

X<1.1)=CMPLX(YC.YS) 

X ( 1 ,NN*?-T ) = X ( 1 « I ) 

continuf 

PFTUPN 

END 


S I ON  AF.XO.GXO.GCONST.WCQNST 

X(2.20A8i,OATA(204yj ,TA(204fl) fA(101 ) ,B (1 01 ) .MFOUP.MF 


COPY  AVAILABLE  TO  DDC  DOES  NOT 
PERMIT  FULLY  LEGIBLE  PRODUCTION 
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SUBROUTINE  EFT  ~(  NSTAGEVSIGN) — 

Complex  x ♦ d at anl » w 
double  precision  d.dt 

DOUPLE  PRECISION  AF , X 0 . 6X0 ♦ GCONST . WCONST 


COMMON  Dt.X (2.2043)  .DAJA  ( 20<*B ) ,TA(204fl)  .A  Cl  01)  »P( 
COMMON  OATANL (2043) .DATAX (2043) 

COMMON  AF(IO) , X 0. GX 0. GCONST ♦ WCONST .NFL AG, I EL AG. M. MM 
INTEGER  R 

N=2<m>NSTAGE  " — 

N2=N/2 

ELTN=N 

PH  1 2N= 6 «.?fl 3 1 S53/FLTM  — * * 

DO  3 J- 1 »~NST AGE  

N2J=N/(?**J) 

NP=N2J 

NI=(?«»J)/2  


DO  2 " I = 1 « NI  

IN?J=( 1-1 ) *N2J 

ELIN2J=IN?J  

TFMP=FLIN?J*PHI?N*SIGN 
W = CMPLX (COS  (TEMP)  ,SIN( 
DO  S R = 1 ♦ NR 
ISUP=R» IN2J 
ISUR1=R*IN2J*2 
lSUP?sISUBl  *N2J 
ISUB3=ISUP*N2 
X (2, ISUa) =X ( 1 . ISUP1 ) *W 
X (?. ISUB3) =X  ( 1 , ISUol ) - 

5 CONTINUE 

2 CONTINUE 
DO  * P = 1 ,N 

X Cl  .R)  =X  (2.R)  " 

6 CONTINUE  

3 CONTINUE 

IF  (SIGN.GT.O.)  pETURN 
^ DO  4 R= 1 ♦ N 
— X (?.R>  =x  (1  .R)/FLTN 

4 CONTINUE 
RETURN 

end  


— SUBROUTINE  TMPROC'  (NF , D , DEL', NST AGE) 

COMPLEX  X. OATANL 

DOUBLE  PRECISION  D.DT  _ 

DOUBLE  PRECISION  AF .XO.GXO. GCONST. WCONST 

COMMON  DT.X (2.2043)  .DATA  ( 2048 ) .TA(20<f3)  , A ( 1 0 1 ) * B C 1 0 1 ) .MFOURtMF 
COMMON  OATANL (2046) .OATAX C204S)  * 

COMMON  AF  (10)  .XO.GXO.GCQNST.’-.'CONST.NFLAG.  IFLAG.M.MM 
LSTAGE=NSTAGE*1 
“ NN=2**LST age 

CALL  SAMPLE  (NF**1 
CALL  FFT  ( LSTAGE 
' DO  2 K=1 » NN 

X ( 1 *K ) =X (2.K) *DATANL (K ) 

2 CONTINUE 

• CALL  FFT  ( LSTAGE.  M.) 

RETURN  ..  

'END  


D.DEL.NSTAGE) 

1.) 


if  it 


I 


FUNCTION  PSI  (D, DEL, SIGN, T) 

COMPLEX  X. DATANL  . 

DOUPLE  PRECISION  D.DT  ' [ 

DOUPLE  PRECISION  AF.x0.GXO. GCONST. wCONST  „ 

COMMON  nT.X<?.?048).DATA(204B).TAl2048),A(101)*B(101)  .MFOUR.MF 
COMMON  DATANL  (204B).i)ATAx(20aA) 


AF (10)  • X 0 . GXO  * GCONST , W CONST ♦ NFL AG » I FLAG, M, MM 
ON  CX(IOI) 


DOUPLE  PRECISION  AF.x0.GXO.Gl 
COMMON  OT,X (2.2048) .DATA (2 
COMMON  DATANL  (?04B)  .L)ATAX  (20 
COMMON  AF  MO)  .XO, GXO. GCONST,' 
DIMENSION  CX(lOl) 

PT=3. 1415927 
Y=2.*Pl*T/D 
C X 1 = COS (Y) 

CX  ( 1 ) =CX  1 

CX  (2) =2 • *C X 1 *CX 1 - 1 . 

DO  2 L=3.MFOUR 

CX (L) =2. *CX CL- 1 ) *CX1-CX (L-2) 

CONTINUE  ' 

z=o. 

DO  3 L=l.MFOUP 
Z=Z-B(L) *CX (L ) /FLOAT (L) 
CONTINUE  *'  ‘ " " 

Z = D«Z 

E=PI*T*T/D 

z=z*e  .. 

E=2.*PI*DFL*T 

G=SIGN«Z 
PSI =G*E 

RETURN * : 

‘END  


SUBROUTINE  TSOPLT 


i X M T N ) 


X MAX=OAT A ( I ) 
XMIN=DATA ( I ) 


DOUPLE  PRECISION  D,DT  7" 1ZI  TIT “771I  T 

‘ DOUPLF  PRECISION  AF • XQ • GXO .GCONST » WCONST  ""  

COMMON  OT.X (2.2048) .DATA (2048) .TA(2048) « A ( 1 0 1 ) « 8(101) »MFOUR»MF 
COMMON  DATANL ( PP4B ) , OAT Ax (2048) 

COMMON  AF (10) »XO«GXo»GCONST .WCONST  » NFL AG. I FLAG »M, MM 
DATA  1ST  AP/ [h*/ 

XMAX=DATA(1) 

XM INsQATA ( 1 ) 

...DO  1 1 = 1 .n  

IF  (QATA(I) .GT.XMAX)  XMAX=UATA(I) 

IF  (DATA ( I ) .LT.XMTN)  XM I N=DAT A ( I ) 

CONTINUE 

oiffsxmax-xmin  

5C= 1 00 . /D IFF  — 

WPITE  (3.300)  SC.XMIN.xv.AX 
> FOPMAT  (4X  . 'SC=‘  .E1P.5,  • mIn= • . E 1 2 . 5 . • 
write  <3.ii7) 

' FOPMAt  ( 1H  ,9H  • - JIME  .5X.1BH  AMPLITUDE  * 

00  1 0 1=1 .N 

XA=APS(SC«(DATA(I)-XMIN) ) 

IDA=IFTX  ( X A ♦ 1 .)  

’ WPITE  (3.200)  TA(I) .DATA(I) , (ISTAR,J=1.IDA) 

I FOPMAT  (1H  .F9.3.5X.E12.5»5X. 102A1 ) 

I CONTINUE 

. |ejurn  .... 


mIn=' .E12.5. • MAX= ' »E 12.5, //////) 
18H  AMPLITUDE  * * ,///)-  "• 


COPY  AVAILABLE  TO  DOG  DOES  NOT 
20  PERMIT  FULLY  LEGIBLE  PRODUCTION 


• *«*  4- 


SUBROUTINE  EVLPSI  (NSTAGE » N » NN> D) 

COMPLEX  X.OATAf.'L 

. r DOUBLE  P>PEClSIOflJ  0«0T  _ ~Z,  _L_ 1."'" 

double  pprcisImm  tr, xo.gxo.gconst.wconst 

DOUBLE  PRECISION  THIGH.OUT 

Common  OT.X  (2.?0<.3)  .0ATA<2ft46)  .TA(?r>4B)  *A(101)  .8(101)  ♦ MFOUR. MF 
COMMON  DAT ANL (204*) .DATAX (?04A) 

...  C0MM0N  aF  (10)  .X0,GX0.6COnST»WCONST*NFLAG»  IFLAG.m.MM  ~ 

TH1GH=0/?. 


I c;y  AGF  = MST  A6F  ♦ 1 

Evaluation  of  the  constants. 


WRITE  (3.1)  ■ - — 

1 format  (ihd 

CALL  INTGPN  (1.0. D 00 , TH I GH , 1 ,5 1 2 ♦ 0 . 1 D-05 , OUT ) 

IF  (NFLAG.GT.O)  GO  TO  BOO  . . 

GCONS1  =OUT  ' * 

CALL  INTGPN  (2.0.0  00,0.51)  00 , 1 . 256 . 0 . ID-05 . OUT ) 

IF  (NFLAG.GT.O)  GO  TO  BBC 

' WCONST=OUT  * 

‘CALL  AMPLE"  (THIGH, N.NN)  

IF  (IFLAG.GT.O. OR. NFLAG.GT.O)  RETURN 
CALL  FFT  (LSTAGF.-l.) 

DO  2 I = 1 , MF OUP  — ~ * 

“ B(I)=-?.*AIMAG(X  (2.1*in  

2 CONTINUE 

DO  10  I = 1 » MFOUR « 2 __  _ 

10  CONTINUE1  — — — " 

RETURN  ..... 

<5<?0  WPITE  n.«) 

;^8  FORMAT  ( 3X , * AN  INTEGRATION  IN  SUBROUTINE  EVLPSI  HAS  FAILED*.//) 

" PFTL'PN  ' ‘ 

end 


SUBROUTINE  SELECT  (NF , VARBLE .OUTPUT)  

COMPLEX  X,  DAT  ANL  

DOUBLE  PRECISION  O.DT  " “ 

DOUBLE  PRECISION  AF  . x 0 • GXO ♦ GCONST » WCOMST 

double  precision  vapble. output. sgnl. taper 

COMMON  OT.X (2. 20Ad) .DATA (20 48) , TA (2048 ) , A < 1 0 1 ) . 

COMMON  DATANL  (20A-3)  .DATAX  (20AB)  ' “ 

COMMON  AF  (10)  .XO.GXO.GCONST.WCONST.NFLAG.IFLAG.M. 
GO  TO  HO, 20)  .NF 

.OUTPUT  = SGNL(VAPBLE)  ‘ 

RETURN  ■'  ' ” ‘ 

0UTPUT=TAOER(VAR8LE)  . 

PF  TURN 

END  • . ,r  s ' -.zmz:  ^ 


8 ( lpl_)  . MFOUR. MF 


COPY  AVAILABLE  TO  DDG  DOES  NOT 
PERMIT  FULLY  LEGIBLE  PRODUCTION 


f\. 


WRKPRoHEP  » c » N1 , NHAX  > STOL ) S J — 

DOUBLE  PRECISION  0»DT 

DOUBLE  RRF C IS  I ON  AF . XO . GXO » GCONST * WCONST  

DOURLE  PRECISION  C.E.STOL.S.DFLOAT.S1 «H .HALF * XY . DS. 
HALF] ,C1 

COMMON  PT.X (2.2048) .DATA (2048) ,TA(2048) «A(101) «B( 
COMMON  OATANL (2048) .DATAX (2048) 

COMMON  AF (10) . XO. 0X0. GCONST. WCONST .NFL AG. IFLAG.M, MM 
NFL A6=0 


— si  =o.oo  oo  — 

CALL  SELECT  (NFiC.FC)  

CALL  SELECT  (NF.E.FE) 

1 M)=?*M] 

IF  (Ml.GT.NMAX)  GO  TO  <399 

H=(F-C) /DFLOAT (Ml)  

HOV2=H/?. 

s=o » D no 

C1=OHOV2  ~ ■'  

“ ‘CALL  SELECT  (NF, Cl. HALF) " 

NM 1 =M 1 - 1 
DO  2 1=1 f NM1 

~ “ XY=C*DFLOAT(I)«H  * 

““'  CALL  SELECT  (NF.XY.S?) 

S=S*S2 
C] =XY+HOV? 

- CALL  SELECT  (NF ,Cl .HALF  1 ) ’ ~ 

HALF=HALF .HALF  1 ~ * - 

2 CONTINUE 

S=(H/A.  ) *(FC->4.»HALF*2.»S*FE) 
" OSsS-Sl 

IF  (OARS (OS) .GT.STOL) ' GO"  TO  3 

RETURN 

3 S)  =5 

GO  TO  1 ■ * ~ 

999  NFLAG  = I — -v~ 

RETURN 

EMO  ' - ' - 


PHASE 


(T. ALPHA) 


“ SUBROUTINE  AMPLE  (THIGH, N,NN)  ~~ ~ 

““COMPLEX  y.OATANL  ~ “ — 

DOUBLE  PPECISIOM  O.DT 

OOURLE  PRECISION  AF  ♦ XO  . GXO.GCONST  . WCONST 

DOUBLE  PRECISION  DFLOAT.  T ♦ T 1 . GT INT . GT INT 1 .Y. RL. THIGH  ‘ 

"COMMON  OT.X (2.2048) ,DATA(?0-fl) .TA(2048) • A ( 1 0 1 ) *8(101) .MFOUR.MF 
CONDON  DATANL (2048) .DATAX (2048) 

COMMON  AF (10) .xo*gxo.gconst«wconst*nflag*iflag*m,mm 

X 0 = 0.0  00  “ - — 

“Gy 0=0.0  oo * * 

GTINT  = 0,6  00  

_ ot*2.*thigh/deloat(nn-1)  “•  ~ ~ ~~  

■“tT=6t/2? “ 

LOOP  FOR  the  EVALUATION  OF  THE  DERIV.  OF  the  PHASE  PSI  (T, ALPHA). 

all  1 NT6RN  c 1 *T  . T 1 . 1 .51  ? * 0.10*05. GT.tNTl.T7  “ — 

IF  (NFLAG.GT.O)  GO  TO  990 
GTINT =GT INT^GT INT 1 

y=wconst*gtint/gconst 

CALL  INVERT  12*0. 10-05. 0. 1D-07,500»Y) ~ — - 

Tr  (NFLAG.GT.O)  GO  TO  990 

IFtJIFLAG.GT.Q)  _RETUR-N _ 

- T ] l K D T ‘ ^ -i  "“ * ~ 

l CONTINUE 
DO  ? 1*1. N 

: rl=(dfloat ( i ) r o . 5i *dt/ (2. ♦thigh):::  . : 

P*DATA(N^I)-R 
X (1 .N.| jsCMPLx (R.0.0) 

* x(i.n-im)=cfplx(-r.6.o)  ~ ■ - 

? CONTINUE  

Rt  TURN 

) WPITE  (3.3) 

j format  ox.BiH  an  integration  in  the  subroutine  sample  has  failed 

l .//> 

return 

ENO 


— ■ '■■?!{g?e,F!;,??o{?^rT'-,N'r,xTou,F-TOL,NTOL,ri--: 

DOURLF  PRECISION  D.OT 

DOUBLE  PRECISION  AF  * XO • GXO « GCONST « WCONST 
IT  . DOUBLE  PRECISION  FX0.X1 .6X1 .XTOL.FTOL.Y.OABS.OERIV.OELTAX 

COMMON  DT.X  (2.2046) .OAT A (2046) .TA(?04b) tA(101).B(101> .MFOUP.MF 
COMMON  OATANL (2046) .OATAX (2046) 

COMMON  AF  (10)  .XO. GXO. GCONST. WCONST.NFLAG.IFLAG.M. MM 

00  1 Is  l »MTOL  * 

FX0=GX0-Y 

IF  (OAHS(FxO) .LT.FTOL)  RETURN 

CALL  SELECT  (NF .XO.DERIV)  ~ " 

’ IF  (OERTV.EC.O.I)  00)  60  TO  999  ~~ 

DELTA X=F  XO/DERlV 
X1=X0 

' X0=X0-PELTAX  ' T’ — - 

CALL  INTGRN  (NF  . X l ,/ 0 .1  .256  ♦ 0. 10-05 . GX  1 ) 

IF  (JjFLAG.C-T.O!  GO  TO  989 
exo=6xo*6* 1 

-T  -KntWu? { DELTAX)  -LT • XT0L > - RETupN-  - 

IF  (OAMS(DELTAX) .GT.XTOL)  GO  TO  990 
RETURN 

989  WRITE  (3.2)  " * — 

— 2 FORMAT  (3X.49H  THE  INTEGRATION  IN  SUBROUTINE  INVERT  HAS  FAILED  * 

RETURN 

999  WRITE  (3.3)  "*  “ - ~ ~ *rwr  

- ~3  FORMAT  (3X.45H  -THE  DERIVATIVE  IN'SU3R0UTINE  ‘INVERT 'IS  ZERO  i//) 

IF LAG=? 

990  WRITE  <T.4) 

4 FORMAT  (3X.50H  DELTAX  JN  SUBROUTINE  INVERT  IS  GREATER  THAN  XTOL 

1 . // ) 

IFLAG=1 

RETURN  ■ __ 


pitpY  AVAILABLE  TO  DD6  DOES  NOT 
PERMIT  TOUT  LEfilBLE  PRODUCTION 
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FUNCTION  TAPER  CPHI) 
iX  x.OATAHL 


SSSSti"  PRECISION  D.OT 

ROUBLE  S»tCf?|<.vj  PI2,CX,YfPHI.CXl,ZtOLOG,OEXP/TAPER 

OOliALt  AF • x 0 » ox o • GCONS T * rtcOMST 

COMMON  r»T  vio  ^ ^ > n*-r  * > «• 

COMMON 
COMMON 


V « 1 4 V”  1 — « " ■"  M » a V » J I f J I 

DT.X  (2.2C48)  .DATA (2048) .TA( 2048)  .A  (101)  *B< 
OATANL  (?04A)  .OATA^  (204=.) 

AF (10) .XO.GXO*GCONST.WCONST»NFLAG,IFLAG*M.MM 

/)Ki  r w # i a t ***  * — * — - ■ * • 


»ocos  : 

101 ) .MFOUR.MF 


"OTMFNSl^N  CX (10) 

P12*6.?8318b307179 

V*PI?<*Pni 

CXlsOCOS(Y) 

- — C»  1 1 ) sCx  l 

CX (?)*?. *CXJ*CX1-1. 

00  ? L*3,S 

CX(L>*?.*CX(L-l)*CXl-CX(L-2) 

2 CONTINUE  "•  

7*1.0  00 
00  3 L*1.S 


(L)  «£X  (L) 


-ONTINUF 

TAPE»=2 

PFTUPN 

END 


function  sgnl  'OT 

COMPLFX  X.OATANL 

DOUPLF  precision  d»dt 

SSSti  J888 ti'oN  -const.  ==■.=:• - 

8«  OA ?I*i <tIif8t!f.i?S5^ii2s?,’TA'so*a>  .«fou«,hf 

COMMON  AF  U 0 ) j X 0,  GX  0 , GCONST  , WCONST  a NFLAG,  I FLAG,  MuMMl'H  'ZZZZ  'JS 
PFTUPN 

ENO  • 


®®#SRS Sff 
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The  purpose  of  this  work  was  to  Identify  and  study  pulse  compression  techniques 
useful  In  LPI  (l.e..  Low  Probability  of  Intercept)  applications.  More  specifically, 
the  technique  proposed  should  be  capable  of  compression  ratios  of  from  128  to  1024 
or  higher,  and  Integrated  sldelobes  of  -35  to  -30  db  throughout  the  compression 
range;  furthermore.  It  should  be  easy  to  mechanize,  and  should  result  In  minimal 
S/N  loss.  It  was  concluded  that  nonlinear  FM  Is  the  pulse-compression  scheme 
capable  of  satisfying  all  these  requirements,  and  existing  technologies  make  Its 
mechanization  feasible. 
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